Finite frequency magnetoelectric response of three dimensional Topological Insulators 
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Topological insulators with a time reversal symmetry breaking perturbation near the surface 
present a magnetoelectric response that is quantized when the frequency of the probing fields is 
much smaller than the surface gap induced by the perturbation. In this work we describe the 
intrinsic finite frequency magnetoelectric response of topological insulators for frequencies of the 
order and larger than the surface gap, including the experimentally relevant case where the system 
is metallic. This response affects physical observable quantities and will give rise to new finite 
frequency phenomena of intrinsic topological origin. 



I. INTRODUCTION 



It is well understood^ that under a time dependent ex- 
ternal electromagnetic perturbation any given material 
will develop a time dependent response, whenever the 
characteristic frequencies of the perturbation arc larger 
than the characteristic frequencies that induce a polar- 
ization or a magnetization in the material. The linear 
response of conventional dielectrics is characterized by 
a dielectric function e(u>) and a magnetic permeability 
/u(w). There is in fact a wider class of materials, known as 
magnetoelectric materials^, which effectively introduce 
other response functions that couple electric and mag- 
netic fields and which in general can also depend on the 
frequency. Perhaps the most striking case of the latter 
are the recently discovered three-dimensional topological 
insulators (TI)&^. They have been predicted to host a 
quantized magnetoelectric term in the action, topologi- 
cal in origin^ of the form S = J dtd 3 x(a/4ir 2 )9~E ■ B, 



where a — e 2 /He and 



only physically relevant 



when time reversal symmetry is broken which implies 
that the surface states that characterize these materials 
are gapped. Despite the fact that this term remains ex- 
perimentally elusive, there has been much ongoing work 
on its consequences. It has has been predicted to give 
rise to a plethora of phenomena including the Kerr and 
Faraday rotation of light determined by the fine struc- 
ture constant^ and a repulsive Casimir effect^iiS, where 
the region of repulsion is determined by 9. Physically, 
these approaches are valid only when the frequencies of 
the relevant fields are much smaller than the surface gap 
m, and the topological magnetoelectric term 6 is inde- 
pendent of the frequency and quantized. In general, the 
magnetoelectric response will depend on the frequency, 
and permeate into physical observablcs, just as the dielec- 
tric function or the magnetic permeability do, modifying 
all the described phenomena related to this topological 
term. 

In this work, we derive the finite frequency magnetoelec- 
tric response of a model hamiltonian which captures the 
basic features of a three dimensional TI. To do so, we 
will generalize the method introduced in Ref. H to finite 
frequency, relating the response of TIs to that of an ef- 



fective model with an extra dimension that behaves as a 
higher dimensional analogue of the Quantum Hall Effect 
(QHE) 1 ^. This approach has been shown to be helpful to 
understand the topological origin of this response, and it 
is also an efficient computational tool in practice. It has 
proven useful to predict the magnetoelectric response in 
a related physical situation^— where the bulk of the 
TI is assumed to be doped. In this particular case, the 
magnetoelectric response is not quantized if the chemical 
potential is outside the band gap, a behaviour that may 
be interpreted as arising from the corresponding anoma- 
lous QHE analogue in five dimensions. 
Here we extend these analysis to a more general and po- 
tentially relevant experimental situation where both the 
frequency and the chemical potential are kept finite, a 
case that can also be understood as descending from a 
five dimensional finite frequency QHE at finite chemi- 
cal potential. As a consistency check we will show how 
known results are recovered in the appropriate limits, 
giving further physical insight into them. 



II. THE MODEL 

Consider the lattice hamiltonian introduced in Ref. [H, 
which captures the low energy description of a generic TI, 
for instance E^Se ^, 13115 - This hamiltonian can be written 
as H = Hq + Hm with 



h.c. - 3t4r°c x , (1) 



(2) 



with x running through all unit cells, s = 1,2,3, and s 
is the lattice vector in the s direction. are defined as 
the set of 4 x 4 matrices that satisfy {r^r^} = 25^, 
with n, v = 0, 1, 2, 3,4 including an extra T4 that will be 
used shortly. In momentum space this can be written as 

H (k) =(tJ2 C0S (M +M-3t)r + ^ sin(/c s )r s 



s=l 



(3) 



2 



These models can be thought of as lattice models 
that host an odd number of low energy massive Dirac 
fermionsAS. For example, for \M\ < t there is a single 
Dirac fcrmion at k = of gap M. 

The presence of a boundary in the hamiltonian can be 
modeled by making the gap position dependent, promot- 
ing © to 



H M = M cos 1 



x)4 r oc x 



(4) 



This accounts for the band inversion by setting 6{— oo) = 
7r in the bulk of the TI and 9{oo) = 0. The specific de- 
pendence of 9 on x will not be needed for our purposes, 
only its asymptotic values. Both experimentally^ and 
from ab initio calculations^ the bulk band gap is well 
approximated by M = 0.3 eV for Bi2Se3. A time rever- 
sal symmetry breaking perturbation may generically be 
included as^ 



H m = 



to sin #(x)c£.r4C x 



(5) 



which is localized at the boundary and opens a surface 
gap to. This surface gap can arise from doping the TI 
with magnetic impurities, and has been measured to be 
to ~ 50 meV-i^. 



III. FINITE FREQUENCY 
ELECTROMAGNETIC RESPONSE OF A 
TOPOLOGICAL INSULATOR 

To obtain the finite frequency response of a TI system 
to electromagnetic fields in the presence of 6{x), we will 
first generalize the original procedure devised in Ref. [B| to 
finite frequency. In what follows, we compute the current 
response of the system with a generalized Kubo formula 
in a way that the effect of 9{x) is included in a manifestly 
perturbative fashion along the derivation. Our starting 
point is to consider the current density at some particular 
point in space-time Xq: 



SS 



SA^xq) ' 



(6) 



where S is the action functional of the system. For a pro- 
file 9{x) that is smooth over length scales l m = 1/{vftti) 
(this is, |V0| << 1/lm): the current at x n is mainly de- 
termined by 9 around 9(xq) = 9q, because correlation 
functions decay exponentially with l m . We may therefore 
include its effects in perturbation theory in 8^9, which is 
by assumption small. For the calculation of j tJ, (xo), we 
thus approximate^: 



(x) « 9(x ) + di9\ x=Xo (x* - 4) + 



(7) 



in the hamiltonian H = H +HM+H m defined by (HI,® 
and ([5]) respectively. To first order in di9 the mass terms 



read 

Hm + H„ 
+ &0L 



^4(Mcos6» r + TOsin6> r 4 )c x 



EC 



x o) c lc(~ M sin 9 T a + to cos #oFi)Cx- 

(8) 



Note that in this hamiltonian 9q is just a constant pa- 
rameter. 

We can now compute the current response at xq when 
a time dependent uniform electric field is applied to the 
system. This is done by computing the expectation value 
of 7 71 to first order in both <9j(9 and the electromagnetic 
field A p , with a generalized Kubo formula 



f(x ) = S 9 



x.x' 



r(x )J^x)JI(x'))A J (x) (9) 



where s = 1,2,3, repeated indices summation is im- 
plied and 20, 2, x' are full space-time variables. In this 
expression J l {x) are the current operators, and Jg{x) is 
the operator attached to d s 9 in © which defines the fol- 
lowing vertex in momentum space 

J$(k) = (-Msin0 o r o + to cos 9 T 4 )d ks = J g d ks . (10) 

With this, the Fourier transform of the current reads 



f(x ) = d s 9 



d 4 p 

BzWf 



TrJl_„/r,Gk-T)Jl „-./iGkJedk 3 Gk 



d 4 k 
zWf 
V <- 



(11) 



J 



where the integral spans the entire Brillouin Zone (BZ). 
The electronic Green's function, which depends on a four- 
momentum vector k = (fco,k) is given by G(k,9 ) = 
(fco — H(k, 9q))~ 1 defined through the Fourier trans- 
formed hamiltonian 

H(k,9 ) = H(k) + (Mcos9 T +msm9 T i ), (12) 

that depends parametrically on 9o, the bulk mass M and 
the surface mass to. The current vertices are defined as 
J* (A;) = 9 %lf o) , J° = I 4 x4 with i = 1,2,3. In the 
derivation we have omitted terms arising from vertices 
with higher derivatives of -ff(k, 9q) with respect to At2£, 
that will not contribute to the magnetoelectric response. 
Before we proceed further, it is worth noting that this 
equation may also be written as 



f(x ) = d s 9\ y 



where 



BZ 



d 4 p 
(27T) 4 



'd„ 



^(P,q,9 ) 



Al 



BZ 



d 4 k 



r ^ T Jk-(p+q)/2Gk-pJi_ p / 2 GkJeGk-q + 



(13) 



(14) 

-p + q 
-> 3 



3 



and again the identity holds disregarding higher deriva- 
tives of H(k, 9 ). II^ (p, q) can be considered the response 
function to 59 = — 6(x) 



BZ 



d 4 p d 4 q 
(2tt) 4 (2tt) 4 



^ + ^n l i{p, q ,e )Aue q , 



(15) 

which is the generalization to finite frequency and mo- 
menta of Ref. 0. Equation represents an equiva- 
lent statement that features an explicit small parameter 
throughout the derivation. 

Consider now the case where the boundary of the TI is in 
the z direction, so that 0(x) = 9(z). A uniform but time 
dependent electric field Ej in momentum space can be 
written in terms of an external vector potential Ai that 
is constant in space, so that Ai(po,p) = <5(p)A;(po)- 
The total current density in the xy plane, shown to 
be quantized in the DC limit£ii2 is defined as J^ D = 
J dzj i (z) with i = x,y. The finite frequency generaliza- 
tion of this quantity, i.e. the integrated current density, 
thus reads 



J2d(po)= / dz d zo 8(z )d qz U^(p,q,9 ) 



Al 



q=0,p=0 



Po ' 



(16) 

With the change of variables /_ dzod Z(i 9 = J_ dd the 
current is finally 



dO 

Ji D (Po)=2n I —d qz irj(p,q,6) 

Z7T 



g=0,p = 



As in the static case, the parameter 9 can be thought of as 
the fifth coordinate of a 4+1 model described by H(k, 9), 
whose response functions are integrated only over half of 
the BZ, because < 9 < tt. The topological finite fre- 
quency response of a 3D TI is thus intimately related to 
the finite frequency response of a D = 4 + 1 insulator. 
Consider now the experimentally relevant limit where 
m <C M. In this limit the response in (fTT|) can be 
obtained analitically. This is due to the fact that the 
low energy physics of hamiltonian (|T2"j) . considered now 
as 4 + 1 hamiltonian in half of the BZ, is dominated 
by an effective 4+1 Dirac fermion of gap m located at 
(k, 9) = (0, 7r/2), where the 4 + 1 analogue of the Berry 
curvature is largest. Integrals in the five dimensional 
BZ are thus well approximated by a region of momenta 
around (0,7r/2) within some cut-off A, and the effective 
hamiltonian in the vicinity of that point is given by 





FIG. 1. Feynman diagrams corresponding to eq. (|19[) . The 
second diagram is obtained from the first by by i O j and 
p <H> — p + q. The grey dot represents the d 3 8 vertex. 



integral in (fl~7|) is 



nfM=/ i^^).. M j (27f)5 



d 5 k 



Tr 



r i G fc _ p P'GfcMr°G fc _, + V^G k+p - q T i G k MY°G 



k — q 



(19) 



where we have used d9 = dk^/M , and eq. (fT4| with the 
current vertices approximated around (0,7r/2): J 1 = 1+ 
J = Mr . The Green functions in these expressions 
are those of a 4+1 Dirac fermion, obtained as G(k) = 
(fco — if (k)) -1 from eq. (fT8)) . where now k has four com- 
ponents. Consequently, the function rdf corresponds 
to the Dirac fermion triangle diagrams shown in Fig. [U 
that can be computed analytically with standard meth- 
ods that are detailed in appendix A (see also Ref. l2ll ). 
and can be considered as the optical response of the five 
dimensional analogue of the QHE. The magnetoelectric 
response is given by the antisymmetric part of the dia- 
gram and thus the total current finally reads 



Ji D (po) = 27rIl6(po,A*)e <3 '^(Po) 
= a(p Q ,fj,)e tJ Ej(p ), 



where the function 



n 5 ( P ,<7) 



2 Po 



q= 0,p=0 



(20) 



(21) 



ipoA J p . 



and we have used that E ] p 

Equations (|20|) and (|2~T|) define the finite frequency re- 
sponse of a TI, which we proceed to evaluate in the next 
section for different experimentally relevant scenarios. 



H(k,9) «/f (0,7r/2) 

=Y l k % + M9T° 



^fci-^l(0,7r/2) k» 

mr 5 



+ d e H \(o,7r/2) 



(18) 



with 9 = 7r/2 + 9. We can identify this hamiltonian as 
that of a 4+1 Dirac fermion where k^ = MO. The cut-off 
for this model is of order A ss M, and for uj « M may 
be taken to infinity. Within this approximation, the 9 



IV. RESULTS 

In this section we compute the function a(p ,fi) de- 
fined above which determines the response of the TI to 
an external electromagnetic field of finite frequency. 
Firstly, it is possible to restore the dependence on the 
chemical potential since nothing in the above argument 



4 




FIG. 2. Real (top) and imaginary (bottom) parts of <r(po, A 1 — 
0) given by (|22[) as a function of po in units of the surface 
gap m. The quantization is lost at higher frequency and a 
logarithmic singularity appears at po = 2m 



depends on whether or not the chemical potential is fi- 
nite as long as \i « M . For a massive Dirac fcrmion at 
zero chemical potential, the function n 5 (p ,/i = 0) can 
be analytically computed. The final analytical expression 
to which one arrives depends only on the surface gap m 
and is given by (we refer the reader to appendix A for 
details) 



1 e 2 to 

a(p ,fi = 0) = --r— log 

2 h po 



2m + p 



2m — po 



(22) 



This function, plotted in Fig. [21 governs the finite fre- 
quency response of a TI and is one of the central results 
of this work. The quantization of the real part at low 
frequencies is broken down at finite frequency giving rise 
to a logarithmic divergence at p = 2m. Therefore, close 
to this range of frequencies, the 9 term will dominate 
the electromagnetic response of the TI. This is particu- 
larly important for Casimir type experiments 9 -^, where 
the interplay between the optical properties of ordinary 
and topological response determines not only the sign of 
the force, but also at what distance does the crossover 
between attractive and repulsive behavior happens. The 
precise way this response alters the Casimir force is an 
interesting issue on its own, and it is left for a future 
study. 

It is important to note that the analytic result (|2"2")l co- 
incides exactly with the optical Hall conductivity of a 
massive D = 2 + 1 dimensional Dirac fermion. The DC 
response of a single Dirac fcrmion is quantized to e 2 /2/i, 
which is consistent with the fact that in the lattice model 
the integrals span only half of the BZ. We thus recover 
the well known result that the boundary of a TI with bro- 
ken time reversal symmetry hosts a half-integer quantum 
Hall effect. 

Being precise, this result should not be interpreted as if 
there is a massive D = 2 + 1 Dirac fermion somewhere in 
the system. Instead, these results imply that the three- 
dimensional optical response of a TI is characterized by 
spatial average in the z direction of all the a xy (z) Hall 



conductivities that occur wherever there is a non zero 
gradient of d z 9(z). This situation is relevant for the re- 
cent experiments described i n 18 ' 19 where TI are doped 
with magnetic impurities that break time reversal sym- 
metry. It is remarkable nevertheless that a full D = 3 + 1 
calculation reduces to a D = 2 + 1 result. As will be 
shown immediately below, this statement does not hold 
for the case of finite chemical potential. 
To do so, one should compute H 5 (p ,^). This can be 
exactly evaluated for some cases which we proceed to de- 
scribe (technical details are left for appendix B). 
One of them is the DC response at finite /i. This limit was 
discussed earlier in KeiJ^— having obvious interest on 
its own since TI appear naturally doped in experiments. 
The numerical evaluation for Hs^po — > 0,/x) in the DC 
limit is exactly given for our model by the expression: 



a(p 0, n) 



sgn(m) 
3m 



IMI' 



if \fj,\ < m 
if > m 



(23) 



The result is shown in Fig. [3ja). The analytical ex- 
pression reveals that although there is a quantized value 
at values of \fi\ < m which also occurs for D = 2 + 1 
fcrmion o 22 ' 23 , already one can notice very important dif- 
ferences with respect to the D = 2 + 1 case. In D = 2 + 1 
it can be show n 22 ! 23 that only a term m/\fx\ arises at fill- 
ings larger than the gap. In the present case however, 
there is a second term which has a different behaviour 
and scales like m 3 /\fj,\ 3 . This term is therefore intrinsi- 
cally related to the D = 3 + 1 nature of the carriers which 
is only fully transparent in the analytic result. The extra 
term turns the kink between the two regimes \fx\ < m 
and |/z| > to smoother, making the curve in Fig. [3ta) 
differcntiable at all fx, in contrast with the 2 + 1 result. 

Finally it is possible to gain analytic insight into the 
regime where both the frequency po and [i are kept fi- 
nite, a situation that can be clearly relevant for experi- 
mentally realistic situations. It is easy to see that when- 
ever < to and po < 2to the result is the same as 
in (|22[) . However there is an experimentally more rele- 
vant situation when > to but still po < 2m. This is 
the case of a doped TI at finite frequency. Evaluating 
Il5(po << 2m, | 1 > to) one obtains the first non trivial 
order in po ■ 



1 e 
1H 



2 r 



3to 



i/'-i 



mp\ 1 
6p 



(24) 



which is plotted as a function of the external frequency 
Po for different values of the chemical potential /x in Fig. 
[3fb). In the limit where (j. = to this coincides with the 
expansion of (f2"2")l when po « 2m. In this general case, 
the quantization of the zero frequency value is also absent 
for finite values of [i and po- Thus our results imply 
that the Kerr and Faraday rotation^ will turn not to 
be quantized in units of the fine structure constant if the 
samples are doped and/or if the frequency of the probe 
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U.J 1 l.J 4 Z.J J , 

(a) (b) 

FIG. 3. a(po, as a function of (a) the chemical potential fi 
for zero frequency in units of the surface gap m and (b) as a 
function of the external frequency po for /i/m = 1.0 — 2.0 in 
steps of 0.2 (top to bottom). There is a quantization plateau 
whenever \fj,\ < m and a decay for /i > m given by (|23|l . For 
finite frequencies and whenever > m is satisfied the DC 
value is not quantized and given by (|24[) . 

is of the order of the surface gap, a common situation in 
actual experiments. 

V. DISCUSSION AND CONCLUSIONS 

To summarize, the present findings will inevitably 
permeate into physically observable quantities whenever 
optical probes have a frequency comparable or larger 
than the surface gap. Given the sizes of the gaps, which 
are as large as 0.3 eV for the bulk gap and 50 meV for 
the surface gap 18 i 19 , it should be possible to observe 
these effects with infra-red probes, which are fully 
controllable within current state of the art technology^. 
The interplay between both scales can be studied in 



full lattice models and will be the aim of a subsequent 
publication. 

More elaborate scenarios, such as the proposed re- 
pulsive Casimir effect£^, the Topological Kerr and 
Faraday effect^ or even the optical-modulator device 
proposed in2£ should be revisited. These findings can 
be generalized to other classes of topological materials 
such as certain classes of Wcyl semi-metals that host 
a Carroll-Field- Jackiw termor— and also to higher 
dimensional analogues of TI. 

In conclusion we have calculated, the electrodynamic 
response of TI at finite frequencies and finite chem- 
ical potential relating it to the response of a higher 
dimensional analogue of the anomalous QHE. Beyond 
reasonable doubt, these findings will permeate and 
strongly affect physical observables, just as any other 
finite frequency response function. We have shown that 
there is a well defined action in this case and that it is 
possible to define a quantity, the total current, which 
is not sensitive to the particular time reversal breaking 
profile inside the TI but it is still dependent on the 
external frequency characterizing the electromagnetic 
perturbation. These results pave the way to the under- 
standing of topological phenomena at finite frequency, 
which are bound to be relevant in current experimental 
set-ups. 
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Appendix A: Calculation of Uz,(po,{i = 0) 

As described in the main text the main text L\^(j>o, A*) determines the finite frequency response of the TI system. In 
this appendix we provide an alternative derivation starting from the response of a D = 4 + 1 system and give details 
of how to compute it for fi = leaving the fi ^ for the last appendix. 

As shown in Ref. [H the quantized DC magnetoelectric response of TI system can be described as descending from a 
five dimensional analogue of the quantized integer quantum Hall effect (IQHE)ii. Under this perspective the work 
of Refs. [T3 - [l4l for finite chemical potential can be understood as arising from the corresponding anomalous QHE 
analogue in five dimensions. Similarly, it is possible to reinterpret our results presented in the main text as descending 
from a five dimensional finite frequency QHE at finite chemical potential. 

To describe the finite frequency response of the D = 4+ 1 at p, = we couple the model to an external electromagnetic 
field A^. After integrating out fcrmions, we obtain an effective action for the gauge field. This effective action will 
generate an analogue of the QHE described by a Chcrn-Simons like term which in momentum space reads: 

^ = J BZ (0 / Wf n 5 (p,q)e^A^A pq<7 A T , (Al) 

where e^ upa is the Levi-Civita totally antisymmetric tensor and A p is the electromagnetic gauge field. In real space, 
the pp, q v momenta turn into derivatives and one recovers an action of the form Ad Ad A, which is a five dimensional 
analogue of the IQHE action in D = 2 + 1 spacctimc dimensions of the form Ad A. The function ^(p, q) accounts 
for the finite frequency, finite momentum response of the system. It is generated in perturbation theory from the 
Feynman diagrams shown in Fig. [T]and can be regarded as arising from the antisymmetric part of the tensor: 

/d 5 k 
j^Tv [G k - p r»G k T»G k - q TP + G k+p - q T v G k T^G k - q TP] . (A2) 

The electronic Green's function is a function of a five-momentum vector k = (fco,k) given by G(k) = (ko — ff(k)) . 
For our model with one massive Dirac fcrmion in D = 4 + 1 dimensions^ the low energy propagator is of the form: 

G(fc) = ^r^- (A3) 

Using that the T a matrices satisfy 

Tr [r r 6 r c r d r e ] = -4e afccde , (A4) 

one can isolate in (|A2[) the antisymmetric term with five r a matrices to obtain, in terms of the Feynman parameters 



it^b, q) = -me'me^p a q T / f dad^ — * (1 ° f l] (A5) 

{^r Jo (k 2 - m 2 +p 2 af3 + q 2 ja + (p + q) 2 /3j) 



= -p^e^ P(T q T fdadftdry -- |L = t ^ ^ = e^ P<7 q T Il 5 (q,p). (A6) 
o 71 " Jo ym — p ap — q ja — (p + qj P7 

It is not difficult to check that this definition of LXs(p, q) is analogous to that given in the main text. To compute it, 
it is possible to numerically evaluate the integrals on the Feynman parameters and find Il^(q,p). 
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As shown in the main text, it is important to keep in mind that in order to calculate the finite frequency response 
to an external time dependent but spatially uniform electric field, only the external frequency po is kept finite while 
all the rest are sent to zero. The integrals in Feynman parameters are analytic and give the logarithmic dependence 
shown in (|22[) in the main text. 

Consistent with the DC response of a single Dirac fermion^, at po = the response is quantized to e 2 /2h, also in 
agreement with the fact that in the lattice model the integrals span only half of the BZ. The theory recovers the fact 
that at the boundary of a TI with broken time reversal symmetry there is a half-integer quantum Hall effect. 



Appendix B: Finite chemical potential: Us(po,n) 



In this appendix we discuss the details of the computation of the response at finite frequency po and chemical 
potential /i. The integral to be computed in this case is defined in (|A5[) with the replacement fco — > fco + M' 



n 5 fo<z) 



-16ie to 



d 5 k 
(2^)5 



dadfid'y 



5(1 



a 



((ko + /-t) 2 k 2 — m 2 + p 2 a(3 + q 2 "fa + (p + q) 2 f3^Y 



(Bl) 



Following the arguments in the main text the relevant case is where all external momenta are zero except po- The 
integral in ko has two third order poles at k^ . The position of the poles in the complex plane is determined by the 
relative magnitude of to 2 , /it, a, po and k 2 . 

Consider the simple case where po —> 0. Following the procedure in Ref. [23| in this case there is a pole which is 
always has a negative imaginary part, no matter what value of k it has. The other pole however depending on k will 
change semi planes and so for certain values of k the integral should be split into two. At this point it is possible to 
identify several cases: 

\/j,\ < to: In this case there are always both poles in different semiplanes and the integral is proportional to sign(m). 

\fi\ > to: In this case it is necessary split the integral on k into two parts. One from to k* and the other one from 
k* to oo where k* is the value of k at which the pole changes semiplane, namely \/ [i 2 — to 2 . The first integral gives 
zero since both poles are on the same side. The second one is: 



n 5 (M) 



16e 2 m 



dko 
27 



k 



dk 



-2tt' 



1 



((k + A*) 2 - k 2 - m 2 )- 



(B2) 



We calculate first the ko integral with the residue theorem. It is a third order pole so closing the contour from above 
and restoring h wc find: 



3 TO 



IMI' 



Since we need cr(/i) = 2^11.5 (//) we finally obtain ([23"1) : 



= 27rn 5 (/i) = - 



3to 



h ' 



(B3) 



(B4) 



which reduces to the familiar i^-sign(?7i) contribution of a 2 + 1 massive Dirac fermion when m = [i, and has an 



extra term 



compared to the D = 2 + 1 casi 



,2:s 



For finite po and \fi\ > m one can generalize the same arguments and find that for po << 2to we have 



a(p < 2m, fx) 



3to 



mpQ 



(B5) 



r 



